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Abstract 

We employ nonperturbative flow equations to compute the equation of state for two 
flavor QCD within an effective quark meson model. This yields the temperature and 
quark mass dependence of quantities like the chiral condensate or the pion mass. A 
precision estimate of the universal critical equation of state for the three-dimensional 
0(4) Heisenberg model is presented. We explicitly connect the 0(4) universal behavior 
near the critical temperature and zero quark mass with the physics at zero temperature 
and a realistic pion mass. For realistic quark masses the pion correlation length near 
T c turns out to be smaller than its zero temperature value. 
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1 Introduction 



Strong interactions in thermal equilibrium at high temperature T — as realized in early 
stages of the evolution of the Universe — differ in important aspects from the well tested 
vacuum or zero temperature properties. A phase transition at some critical temperature T c 
or a relatively sharp crossover may separate the high and low temperature physics [1]. Many 
experimental activities at heavy ion colliders [2] search for signs of such a transition. It was 
realized early that the transition should be closely related to a qualitative change in the chiral 
condensate according to the general observation that spontaneous symmetry breaking tends 
to be absent in a high temperature situation. A series of stimulating contributions [3, 4, 5] 
pointed out that for sufficiently small up and down quark masses, m u and m d , and for a 
sufficiently large mass of the strange quark, m s , the chiral transition is expected to belong 
to the universality class of the 0(4) Heisenberg model. This means that near the critical 
temperature only the pions and the sigma particle play a role for the behavior of condensates 
and long distance correlation functions. It was suggested [4, 5] that a large correlation length 
may be responsible for important fluctuations or lead to a disoriented chiral condensate [6]. 
This was even related [4, 5] to the spectacular "Centauro events" [7] observed in cosmic rays. 
The question how small m u and m d would have to be in order to see a large correlation length 
near T c and if this scenario could be realized for realistic values of the current quark masses 
remained, however, unanswered. The reason was the missing link between the universal 
behavior near T c and zero current quark mass on one hand and the known physical properties 
at T = for realistic quark masses on the other hand. 

It is the purpose of the present paper to provide this link. We present here the equation 
of state for two flavor QCD within an effective quark meson model. The equation of state 
expresses the chiral condensate (ijjifj as a function of temperature and the average current 
quark mass rh = (m u + m d )/2. This connects explicitly the universal critical behavior for 
T — > T c and rh — > with the temperature dependence for a realistic value m p h ys - Since 
our discussion covers the whole temperature range < T^1.7T C we can fix m p h ys such 
that the (zero temperature) pion mass is m n = 135 MeV. The condensate (ipip^ plays here 
the role of an order parameter. Its precise definition will be given in section 2. Figure 1 
shows our results for (ipip^ (T, rh): Curve (a) gives the temperature dependence of (V^) 
in the chiral limit rh — 0. Here the lower curve is the full result for arbitrary T whereas 
the upper curve corresponds to the universal scaling form of the equation of state for the 
0(4) Heisenberg model. We see perfect agreement of both curves for T sufficiently close 
to T c = 100.7 MeV. This demonstrates the capability of our method to cover the critical 
behavior and, in particular, to reproduce the critical exponents of the 0(4)-model. We have 
determined (cf. section 5) the universal critical equation of state as well as the non-universal 
amplitudes. This provides the full functional dependence of {T, 171 ) f° r small T — T c 

and rh. The curves (b), (c) and (d) are for non-vanishing values of the average current 
quark mass rh. Curve (c) corresponds to m phys or, equivalently, m n (T = 0) = 135 MeV. 
One observes a crossover in the range T = (1.2 — 1.5)T C . The 0(4) universal equation 
of state (upper curve) gives a reasonable approximation in this temperature range. The 
transition turns out to be much less dramatic than for rh — 0. We have also plotted in 
curve (b) the results for comparably small quark masses ~ 1 MeV, i.e. rh = m phys /10, for 
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Figure 1: The plot shows the chiral condensate (ipip) as a function of temperature T. Lines (a), (6), (c), 
(e?) correspond at zero temperature to = 0, 45 MeV, 135 MeV, 230 MeV, respectively. For each pair of 
curves the lower one represents the full T-depcndcncc of (ipip) whereas the upper one shows for comparison 
the universal scaling form of the equation of state for the 0(4) Hcisenbcrg model. The critical temperature 
for zero quark mass is T c = 100.7 MeV. The chiral condensate is normalized at a scale /c$ ~ 620 MeV. 

which the T = value of m n equals 45 MeV. The crossover is considerably sharper but 
a substantial deviation from the chiral limit remains even for such small values of rh. In 
order to facilitate comparison with lattice simulations which are typically performed for 
larger values of we also present results for m„-(T = 0) = 230 MeV in curve (d). One 
may define a "pseudocritical temperature" T pc associated to the smooth crossover as the 
inflection point of (^>^ (T) as usually done in lattice simulations. Our results for this 

definition of T pc are denoted by T^) and are presented in table 1 for the four different 
values of m or, equivalently, m n (T = 0). The value for the pseudocritical temperature for 
= 230 MeV compares well with the lattice results for two flavor QCD (cf. section 5). One 
should mention, though, that a determination of T pc according to this definition is subject 
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Table 1: The table shows the critical and "pseudocritical" temperatures for various values of the zero 
temperature pion mass. Here is defined as the inflection point of ^tpip) (T) whereas T$ is the location 
of the maximum of the sigma correlation length (see section 4). 
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Figure 2: The plot shows as a function of temperature T for three different values of the average light 
current quark mass rh. The solid line corresponds to the realistic value rh = m p h ys whereas the dotted line 
represents the situation without explicit chiral symmetry breaking, i.e., rh = 0. The intermediate, dashed 
line assumes rh = m p h ys /10. 

to sizeable numerical uncertainties for large pion masses as the curve in figure 1 is almost 
linear around the inflection point for quite a large temperature range. A problematic point 
in lattice simulations is the extrapolation to realistic values of m n or even to the chiral limit. 
Our results may serve here as an analytic guide. The overall picture shows the approximate 
validity of the 0(4) scaling behavior over a large temperature interval in the vicinity of and 
above T c once the (non-universal) amplitudes are properly computed. 

A second important result of our investigations is the temperature dependence of the 
space-like pion correlation length m~ l (T). (We will often call m n (T) the temperature 
dependent pion mass since it coincides with the physical pion mass for T = 0.) The plot for 
m n (T) in figure 2 again shows the second order phase transition in the chiral limit rh — 0. 
For T <T C the pions are massless Goldstone bosons whereas for T > T c they form with the 
sigma particle a degenerate vector of 0(4) with mass increasing as a function of temperature. 
For m = the behavior for small positive T — T c is characterized by the critical exponent 
v, i.e. m n (T) = (^ + )" 1 T C ((T - T c )/T c ) u and we obtain v = 0.787, = 0.270. For rh > 
we find that m^T) remains almost constant for T ^ T c with only a very slight dip for T 
near T c /2. For T > T c the correlation length decreases rapidly and for T>T C the precise 
value of rh becomes irrelevant. We see that the universal critical behavior near T c is quite 
smoothly connected to T = 0. The full functional dependence of m n (T,rh) allows us to 
compute the overall size of the pion correlation length near the critical temperature and we 
find m n (T pc ) ~ 1.7m 7r (0) for the realistic value rh phys . This correlation length is even smaller 
than the vacuum (T = 0) one and gives no indication for strong fluctuations of pions with 
long wavelength. It would be interesting to see if a decrease of the pion correlation length 
at and above T c is experimentally observable. It should be emphasized, however, that a 
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tricritical behavior with a massless excitation remains possible for three flavors. This would 
not be characterized by the universal behavior of the 0(4)-model. We also point out that 
the present investigation for the two flavor case does not take into account a speculative 
"effective restoration" of the axial Ua(X) symmetry at high temperature [3, 8]. We will 
comment on these issues in section 6. 

Our method is based on the effective average action T k [9] which is the generating func- 
tional of the 1PI Green functions in presence of an infrared cutoff k. In a thermal equilib- 
rium context Tfc depends also on temperature and describes a coarse grained free energy as 
a functional of appropriate fields. Here k~ l corresponds to the coarse graining length scale. 
Varying the infrared cutoff k allows us to consider the relevant physics in dependence on 
some momentum-like scale. The results for the order parameter and correlation functions 
presented in this paper are obtained by removing the infrared cutoff (k — ■> 0) in the end. 
For scalar fields $j the /c-dependence of the effective average action is given by an exact 
nonperturbative flow equation [10] 



d_ 
di 



r k m = \Trl [ (T^m + R k )' ld ^ (1.1) 



where t = ln(/c/A) with A an arbitrary momentum scale. Here rj^ denotes the matrix of 
second functional derivatives of T k with respect to the field components: 



and we employ a momentum dependent infrared cutoff 

Z^e-« 2 / k2 
= x _ e - q y k i ( L3 ) 

with Z^ k an appropriate wave function renormalization constant to be specified later. In 
momentum space the trace contains a momentum integration, Tr = / The flow 

(2) 

equation (1.1) closely resembles a one-loop equation: Indeed, replacing T k by the second 

(2) 

functional derivative of the classical action, , it corresponds to the one-loop result for a 
theory where an infrared cutoff 

I J (^d^te)flk(g)¥>i(-g) 

is added to the classical action S c \ [</?]■ This cutoff appears in the inverse "average propaga- 
tor" 

P(q) = q 2 + Z^ k R k (q) = ^L-^ (1.4) 

l-exp{-fj) 

which approaches k 2 for q 2 <ti k 2 . Up to exponentially small corrections the integration 
of the high momentum modes with q 2 ^> k 2 is not affected by the infrared cutoff. The 
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"renormalization group improvement" S£i ~> contains all contributions beyond one- 
loop and makes (1.1) exact. Of course, it also turns the flow equation into a functional 
differential equation which cannot be solved exactly in general. We emphasize that the 
flow equation (1.1) is connected to the Wilsonian renormalization group equation [11]— [15] 
(often also called exact renormalization group equation). Extensions of the flow equations 
to fermions [16, 17] and gauge fields [18]-[24] are available. 

Since in most cases the flow equation can not be solved exactly the capacity to devise 
useful truncations in a nonperturbative context becomes crucial. This requires first of all 
an identification of the degrees of freedom which are most relevant for a given problem. 
In the present paper we concentrate on the chiral aspects of QCD 1 . Spontaneous chiral 
symmetry breaking occurs through the expectation value of a (complex) scalar field $ a b 
which transforms as (N, N) under the chiral flavor group SUl(N) x SUr(N) with N the 
number of light quark flavors. More precisely, the expectation value 

<V b ) = cf <5 a6 (1.5) 

induces for Wq ^ a spontaneous breaking of the chiral group to a vector-like subgroup, 
SUl(N) x SUr(N) — ► SUl+r(N) = SUv(N). In addition, non-vanishing current quark 
masses m u ,md,m s break the chiral group explicitly and also lift the SUv(N) degeneracy of 
the spectrum if they are unequal. The physical degrees of freedom contained in the field 
$ a b are pseudoscalar and scalar mesons which can be understood as quark-ant iquark bound 
states. It is obvious that any analytical description of the chiral transition has to include 
at least part of these (pseudo-) scalar fields as the most relevant degrees of freedom. 

In the present work we use for k smaller than a "compositeness scale" k$ ~ 600 MeV a 
description in terms of <3> a ft and quark degrees of freedom. The quarks acquire a constituent 
quark mass M q through the chiral condensate Zfo which forms in our picture for k x sB — 
400 MeV. This effective quark meson model can be obtained from QCD by "integrating 
out" the gluon degrees of freedom and introducing fields for composite operators [26, 22]. 
This will be explained in more detail in the first part of section 2. In this picture the scale 
fc$ is associated to the scale at which the formation of mesonic bound states can be observed 
in the flow of the effective (momentum dependent) four-quark interaction. We will restrict 
our discussion in this paper to two flavor QCD with equal quark masses m u = = rh. 
Since in this case the scalar triplet a and the pseudoscalar singlet (associated with the rf) 
have typical masses around 2 1 GeV we will neglect them for k < This reduces the scalar 
degrees of freedom of our effective model to a four component vector of 0(4), consisting of 
the three pions and the "sigma resonance" . 

We imagine that all other degrees of freedom besides the quarks ip and the scalars $ are 
integrated out. This is reflected in the precise form of the effective average action [ip, $] 
at the scale k$ which serves as an initial value for the solution of the flow equation. The flow 

■"■For a study of chiral symmetry breaking in QED using related exact renormalization group techniques 
see [25]. 

2 More precisely, because of the anomalous Ua{1) breaking in QCD these mesons are significantly heavier 
than the remaining degrees of freedom in the range of scales k where the dynamics of the model is strongly 
influenced by mesonic fluctuations. The situation becomes more involved if the model is considered at high 
temperature which is discussed in section 6. 
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of for k < k<$> is then entirely due to the quark and meson fluctuations which are 

not yet included in [ip, $]. Obviously, the initial value may be a quite complicated 
functional of ip and $ containing, in particular, important non-local behavior. We will 
nevertheless use a rather simple truncation in terms of standard kinetic terms and a most 
general form of the scalar potential U k) i.e. 3 



Here T k is invariant under the chiral flavor symmetry SUl(2) x SUr(2) and the only explicit 
symmetry breaking arises through the source term j ~ rh. We will consider the flow of 
the most general form of Uk consistent with the symmetries (without any restriction to 
a polynomial form as typically used in a perturbative context). On the other hand, our 
approximations for the kinetic terms are rather crude and parameterized by only two running 
wave function renormalization constants, and Z^^- The same holds for the effective 
Yukawa coupling h k . The main approximations in this work concern 

(i) the simple form of the derivative terms and the Yukawa coupling, in particular, the 

neglect of higher derivative terms (and terms with two derivatives and higher powers 
of $). This is partly motivated by the observation that at the scale k$ and for small 
temperatures the possible strong non-localities related to confinement affect most 
likely only the quarks in a momentum range q 2 ^ (300 MeV) 2 . Details of the quark 
propagator and interactions in this momentum range are not very important in our 
context (see section 2). 

(ii) the neglect of interactions involving more than two quark fields. This is motivated by 

the fact that the dominant multi-quark interactions are already incorporated in the 
mesonic description. Six-quark interactions beyond those contained effectively in Uk 
could be related to baryons and play probably only a minor role for the meson physics 
considered here. 

We will choose a normalization of if},<& such that Z^^ — = 1. We therefore need 
as initial values at the scale k$ the scalar wave function renormalization Z$ t k* an d the 
shape of the potential U^. We will make here the important assumption that Z$^ is 
small at the compositeness scale k$ (similarly to what is usually assumed in Nambu-Jona- 
Lasinio-like models). This results in a large value of the renormalized Yukawa coupling 
hk = Z^ k Z^ k hk- A large value of is phenomeno logically suggested by the comparably 
large value of the constituent quark mass M q . The latter is related to the value of the 
Yukawa coupling for k — > and the pion decay constant f w = 92.4 MeV by M q = hf w /2 

3 Our Euclidean conventions (hk is real) are specified in refs. [27, 16]. 





(1.6) 
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(with h = hk=o), and M q ~ 300 MeV implies h 2 /An ~ 3.4. For increasing k the value of 
the Yukawa coupling grows rapidly for k ^ M q . Our assumption of a large initial value 
for is therefore equivalent to the assumption that the truncation (1.6) can be used 
up to the vicinity of the Landau pole of h^. The existence of a strong Yukawa coupling 
enhances the predictive power of our approach considerably. It implies a fast approach of 
the running couplings to partial infrared fixed points [27]. In consequence, the detailed form 
of Uk^, becomes unimportant, except for the value of one relevant parameter corresponding 
to the scalar mass term ra| . In this paper we fix such that f n = 92.4 MeV for 
m w = 135 MeV. The possibility of such a choice is highly non-trivial since f n can actually 
be predicted [27] in our setting within a relatively narrow range. The value f n = 92.4 MeV 
(for m n = 135 MeV) sets our unit of mass for two flavor QCD which is, of course, not directly 
accessible by observation. Besides Tri\ (or /„.) the other input parameter used in this work 
is the constituent quark mass M q which determines the scale at which becomes very 
large. We consider a range 300 MeV ^ M q ^ 350 MeV and find a rather weak dependence 
of our results on the precise value of M q . We also observe that the limit h ks> — > oo can be 
considered as the lowest order of a systematic expansion in which is obviously highly 
nonperturbative. 

A generalization of our method to the realistic case of three light flavors is possible and 
work in this direction is in progress. For the time being we expect that many features found 
for iV = 2 will carry over to the realistic case, especially the critical behavior for T — > T c and 
rh — > (for fixed m s ^ 0). Nevertheless, some quantities like (v/iip^ (T — 0), the difference 
between f n for realistic quark masses and rh = or the mass of the sigma resonance at 
T = may be modified. This will also affect the non-universal amplitudes in the critical 
equation of state and, in particular, the value of T c . In the picture of the two flavor quark 
meson model these changes occur through an effective temperature dependence of the initial 
values of couplings at the scale This effect, which is due to the temperature dependence 
of effects from fluctuations not considered in the present work is discussed briefly in section 
6. It remains perfectly conceivable that this additional temperature dependence may result 
in a first order phase transition or a tricritical behavior for realistic values of rh for the three 
flavor case. Details will depend on the strange quark mass. We observe, however, that the 
temperature dependence in the limit rh — ■> involves for T <T C only information from the 
running of couplings in the range k ^ 300 MeV. (The running for k ^ 3T effectively drops 
out in the comparison between the thermal equilibrium results and those for T = 0.) In 
this range of temperatures our model should be quite reliable. 

Finally, we mention that we have concentrated here only on the ^-dependent part of the 
effective action which is related to chiral symmetry breaking. The ^-independent part of 
the free energy also depends on T and only part of this temperature dependence is induced 
by the scalar and quark fluctuations considered in the present paper. Most likely, the gluon 
degrees of freedom cannot be neglected for this purpose. This is the reason why we do not 
give results for "overall quantities" like energy density or pressure as a function of T. 

This paper is organized as follows: In section 2 we review the linear quark meson model 
at vanishing temperature. We begin with an overview of the different scales appearing in 
strong interaction physics. Subsequently, the flow equations for the linear quark meson 
model are introduced and their approximate partial fixed point behavior is discussed in 



7 



detail leading to a "prediction" of the chiral condensate In section 3 the exact 

renormalization group formulation of field theories in thermal equilibrium is given. It is 
demonstrated how mass threshold functions in the flow equations smoothly decouple all 
massive Matsubara modes as the temperature increases, therefore leading to a "dimensional 
reduction" of the model. Section 4 contains our results for the linear quark meson model 
at non-vanishing temperature. Here we discuss the T-dependences of the parameters and 
physical observables of the linear quark meson model in detail for a temperature range 
< T ^ 170 MeV including the (pseudo) critical temperature T c of the chiral transition. 
The critical behavior of the model near T c and m — 0, where m denotes the light average 
current quark mass, is carefully analyzed in section 5. There we present the universal scaling 
form of the equation of state including a semi-analytical fit for the corresponding scaling 
function. Also the universal critical exponents and amplitude ratios are given there. The 
effects of additional degrees of freedom of strong interaction physics not included in the linear 
0(4)-symmetric quark meson model are addressed in section 6. Here we also comment on 
differences between the linear quark meson model and chiral perturbation theory. Some 
technical details concerning the quark mass term and the definition of threshold functions 
at vanishing and non-vanishing temperature are presented in three appendices. 

2 The quark meson model at T = 

Before discussing the finite temperature behavior of strong interaction physics we will re- 
view some of its zero temperature features. This will be done within the framework of a 
linear quark meson model as an effective description for QCD for scales below the mesonic 
compositeness scale of approximately k$ ~ 600 MeV. Relating this model to QCD in a 
semi-quantitative way in subsection 2.1 will allow us to gain some information on the initial 
value for the effective average action at the compositeness scale k®. We emphasize, however, 
that the quantitative aspects of the derivation of the effective quark meson model from QCD 
will not be relevant for our practical calculations in the mesonic sector. This is related to 
the "infrared stability" for large Yukawa coupling as discussed in the introduction and 
which will be made quantitative in subsection 2.2. 

2.1 A short (scale) history of QCD 

For an evaluation of the trace on the right hand side of the flow equation (1.1) only a 
small momentum range q 2 ~ k 2 contributes substantially. One therefore only needs to take 
into account those fluctuations which are important in this momentum interval. Here we 
are interested in the description of chiral symmetry breaking. The relevant fluctuations in 
relation to this phenomenon may change with the scale k and we begin by summarizing 
the qualitatively different scale intervals which appear for meson physics in QCD. Some of 
this will be explained in more detail in the remainder of this section whereas other aspects 
are well known. Details of this discussion may also be found in refs. [26, 27, 28]. We will 
distinguish five qualitatively different ranges of scales: 
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1. At sufficiently high momentum scales, say, 

k £ k p ~ 1.5 GeV 

the relevant degrees of freedom of strong interactions are quarks and gluons and their 
dynamics is well described by perturbative QCD. 

2. For decreasing momentum scales in the range 

k$ ~ 600 MeV £ A; £ jfep ~ 1.5 GeV 

the dynamical degrees of freedom are still quarks and gluons. Yet, as k is lowered part 
of their dynamics becomes dominated by effective non-local four quark interactions 
which cannot be fully accessed perturbatively. 

3. At still lower scales this situation changes dramatically. Quarks and gluons are 
supplemented by mesonic bound states as additional degrees of freedom which are 
formed at a scale k<$> ~ 600 MeV. We emphasize that kq> is well separated from 
Aq CD ~ 200 MeV where confinement sets in and from the constituent masses of the 
quarks M q ~ (300 — 350) MeV. This implies that below the compositeness scale /c$ 
there exists a hybrid description in term of quarks and mesons! It is important to note 
that for scales not too much smaller than chiral symmetry remains unbroken. This 
situation holds down to a scale k x sB — 400 MeV at which the scalar meson potential 
develops a non-trivial minimum thus breaking chiral symmetry spontaneously. The 
meson dynamics within the range 

k xSB ^ 400 MeV < k < k$ ~ 600 MeV 

is dominated by light current quarks with a strong Yukawa coupling hl/(Air) a s (k) 
to mesons. We will thus assume that the leading gluon effects are included below k$ 
already in the formation of mesons. Near k x $B also fluctuations of the light scalar 
mesons become important as their initially large renormalized mass approaches zero. 
Other hadronic bound states like vector mesons or baryons should have masses larger 
than those of the lightest scalar mesons, in particular near k x sB, and give therefore 
only subleading contributions to the dynamics. This leads us to a simple linear model 
of quarks and scalar mesons as an effective description of QCD for scales below 

4. As one evolves to scales below k x ss the Yukawa coupling decreases whereas a s in- 
creases. Of course, getting closer to Aqcd it is no longer justified to neglect in the 
quark sector the QCD effects which go beyond the dynamics of the effective quark me- 
son model in our truncation (1.6). On the other hand, the final IR value of the Yukawa 
coupling h is fixed by the typical values of constituent quark masses M q ~ 300 MeV 
to be h 2 /(Air) ~ 3.4. One may therefore speculate that the domination of the Yukawa 
interaction persists even for the interval 

M q ~ 300 MeV £ k £ A; xSB ^ 400 MeV 
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below which the quarks decouple from the evolution of the mesonic degrees of freedom 
altogether. Of course, details of the gluonic interactions are expected to be crucial 
for an understanding of quark and gluon confinement. Strong interaction effects may 
dramatically change the momentum dependence of the quark propagator for k and q 2 
around A QCD . Yet, there is no coupling of the gluons to the color neutral mesons. As 
long as one is only interested in the dynamics of the mesons one is led to expect that 
confinement effects are quantitatively not too important. 

5. Because of the effective decoupling of the quarks and therefore of the whole colored 
sector the details of confinement have only little influence on the mesonic dynamics 
for scales 



Here quarks and gluons disappear effectively from the spectrum and one is left with 
the pions. They are the only particles whose propagation is not suppressed by a large 
mass. For scales below the pion mass the flow of the couplings stops. 

Let us now try to understand these different ranges of scales in more detail. We may start 
at k p — 1.5 GeV where we assume that all gluonic degrees of freedom have been integrated 
out while we have kept an effective infrared cutoff ~ k p in the quark propagators. Details of 
this procedure were outlined in ref . [22] . This results in a non-trivial momentum dependence 
of the quark propagator and effective non-local four and higher quark interactions. Because 
of the infrared cutoff the resulting effective action for the quarks resembles closely the one 
for heavy quarks (at least for Euclidean momenta). The dominant effect is the appearance 
of an effective quark potential similar to the one for the charm quark which describes the 
effective four quark interactions. 

We next have to remove the infrared cutoff for the quarks, k — > 0. This task can be 
carried out by means of the exact flow equation for quarks only, starting at k p with an 
initial value [^J as obtained after integrating out the gluons. A first investigation in 
this direction [26] used a truncation with a chirally invariant four quark interaction whose 
most general momentum dependence was retained. A crucial point is, of course, the initial 
value for this momentum dependence at k p . The ansatz used in ref. [26] is obtained by 
Fierz transforming the heavy quark potential and keeping, for simplicity, only the scalar 
meson channel while neglecting the p-meson and pomeron channels which are also present. 
The effective heavy quark potential was approximated there by a one gluon exchange term 
~ a s {k p ) supplemented by a linearly rising string tension term. This ansatz corresponds 
to the four quark interaction generated by the flavor neutral t-channel one gluon exchange 
depicted in figure 3 with an appropriately modified gluon propagator and quark gluon vertex 
in order to account for the linearly rising part of the potential. 

The evolution equation for the four quark interaction can be derived from the fermionic 
version of eq. (1.1). It is by far not clear that the evolution of the effective four quark vertex 
will lead at lower scales to a momentum dependence representing the (s-channel) exchange 
of colorless mesonic bound states. Yet, at the compositeness scale 



k £ M Q ~ 300 MeV . 



k«, ~ 600 MeV 
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fcp = 1.5GeV A; $ ~ 600MeV 

Figure 3: The left diagram represents the one gluon exchange ^-channel contribution to the four quark 
vertex at the scale k p ~ 1.5 GeV. It is assumed here that the gluon propagator is modified such that it 
accounts for the linearly rising term in the heavy quark potential. The right diagram displays the scalar 
meson s-channel exchange found at the compositeness scale fc$ ~ 600 MeV. 

one finds [26] an approximate Bethe-Salpeter factorization of the four quark amplitude 
with precisely this property. This situation is described by the right Feynman diagram 
in figure 3. In particular, it was possible to extract the amputated Bethe-Salpeter wave 
function as well as the mesonic bound state propagator displaying a pole-like structure in 
the s-channel if it is continued to negative s = (pi +P2) 2 - In the limit where the momentum 
dependence of the Bethe-Salpeter wave function and the bound state propagator is neglected 
the effective action T^ resembles 4 the Nambu-Jona-Lasinio model [29, 30]. It is therefore 
not surprising that our description of the dynamics for k < k$, will parallel certain aspects of 
the investigations of this model, even though we are not bound to the approximations used 
typically in such studies (large- N c expansion, perturbative renormalization group, etc.). 

It is clear that for scales k ^ kq> a description of strong interaction physics in terms 
of quark fields alone would be rather inefficient. Finding physically reasonable truncations 
of the effective average action should be much easier once composite fields for the mesons 
are introduced. The exact renormalization group equation can indeed be supplemented by 
an exact formalism for the introduction of composite field variables or, more generally, a 
change of variables [26]. In the context of QCD this amounts to the replacement of the 
dominant part of the four quark interactions by scalar meson fields with Yukawa couplings 
to the quarks. In turn, this substitutes the effective quark action at the scale k<$> by the 
effective quark meson action given in eq. (1.6) in the introduction 5 . The term in the effec- 

4 Our solution of the flow equation for Tk with Z^^* = (see below) may be considered as a solution of 
the NJL model with a particular form of the ultraviolet cutoff dictated by the shape of Rk(q 2 ) as given in 
cq. (1.3). 

5 We note that no double counting problem arises in this procedure. 
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tive potential which is quadratic in <3>, = rn 2 k<6 tr + . . ., turns out to be positive 
as a consequence of the attractiveness of the four quark interaction inducing it. Its value 
was found for the simple truncations used in ref. [26] to be ~ 120 MeV. The higher 
order terms in cannot be determined in the four quark approximation since they corre- 
spond to terms involving six or more quark fields. (Their values will not be needed for our 
quantitative investigations as is discussed in subsection 2.2). The initial value of the (bare) 
Yukawa coupling corresponds to the amputated Bethe-Salpeter wave function. Neglecting 
its momentum dependence it can be normalized to = 1. Moreover, the quark wave 
function renormalization is normalized to one at the scale k§ for convenience. One 
may add that we have refrained here for simplicity from considering four quark operators 
with vector and pseudo-vector spin structure. Their inclusion is straightforward and would 
lead to vector and pseudo-vector mesons in the effective action. 

In view of the possible large truncation errors made in ref. [26] we will take (2.1) and the 
above value of fn^ only as order of magnitude estimates. Furthermore, we will assume, as 
motivated in the introduction and usually done in large- N c computations within the NJL 
model, that 

Z* M < 1 . (2.2) 

1 /2 l 

As a consequence, the initial value of the renormalized Yukawa coupling — Z^ k ^ 
hk 9 is much larger than one and we will be able to exploit the infrared stable features of the 
flow equations. As a typical coupling we take = 100 in order to simulate the limit — > 
oo. The effective potential £4(4>) must be invariant under the chiral SUl{N)xSUr(N) flavor 
symmetry. In fact, the axial anomaly of QCD breaks the Abelian 6^(1) symmetry. The 
resulting Ua(X) violating multi-quark interactions 6 lead to corresponding Ua(X) violating 
terms in £4($). Accordingly, the most general effective potential U k is a function of the 
N + 1 independent C and P conserving SU L (N) x SU R (N) invariants 

p = tr $ t( I> , 

n ~ tr(V$_I p y , i = 2,...,N, (2.3) 
f = det $ + det $ f . 

For a given initial form of Uk all quantities in our truncation of Tk (1-6) are now fixed and 
we may follow the flow of Tk to k — > 0. In this context it is important that the formalism 
for composite fields [26] also induces an infrared cutoff in the meson propagator. The flow 
equations are therefore exactly of the form (1.1), with quarks and mesons treated on an 
equal footing. At the compositeness scale the quadratic term of Uk 9 = m| Tr + ... is 
positive and the minimum of therefore occurs for $ = 0. Spontaneous chiral symmetry 
breaking is described by a non-vanishing expectation value (<&) in absence of quark masses. 
This follows from the change of the shape of the effective potential Uk as k flows from 
to zero. The large renormalized Yukawa coupling rapidly drives the scalar mass term 
to negative values and leads to a potential minimum away from the origin at some scale 

6 A first attempt for the computation of the anomaly term in the fermionic effective average action can 
be found in ref. [31]. 
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k x SB < kq> such that finally (<£>) = a 7^ for k — > [26, 27]. This concludes our overview of 
the general features of chiral symmetry breaking in the context of flow equations for QCD. 

We will concentrate in this work on the two flavor case (N = 2) and comment on the 
effects of including the strange quark in section 6. Furthermore we will neglect isospin 
violation and therefore consider a singlet source term j proportional to the average light 
current quark mass rh = \{m u + ma). Due to the C/A(l)-anomaiy there is a mass split for 
the mesons described by $. The scalar triplet (a ) and the pseudoscalar singlet (?/) receive 
a large mass whereas the pseudoscalar triplet (ir) and the scalar singlet (a) remains light. 
From the measured values m v r,m ao ~ 1 GeV it is evident that a decoupling of these mesons 
is presumably a very realistic limit 7 . It can be achieved in a chirally invariant way and 
leads to the well known 0(4)-symmetric Gell-Mann-Levy linear sigma model [32] which is, 
however, coupled to quarks now. This is the two flavor linear quark meson model which 
we will study in the remainder of this work. For this model the effective potential U k is a 
function of p only. 

It remains to determine the source j as a function of the average current quark mass 
rh. This is carried out in appendix A and we obtain in our normalization with = 1, 

hk 9 = 1, 

J = 2m^m . (2.4) 

It is remarkable that higher order terms do not influence the relation between j and rh. 
Only the quadratic term m£ enters which is in our scenario the only relevant coupling. 
This feature is an important ingredient for the predictive power of the model as far as the 
absolute size of the current quark mass is concerned. 

The quantities which are directly connected to chiral symmetry breaking depend on the 
fc-dependent expectation value ($) fe = a 0jk as given by 



(P = 2a") = (2.5) 



In terms of the renormalized expectation value 

°o,fc = Z^lao,k (2.6) 
we obtain the following expressions for phenomenological observables from (1.6) for 8 d — 4 

fir,k = 2<T fc , 



1/2 

- m 



2 7 -l/2 ^j> _ 7 -l/2 J 

&0,k 2(T 0)fc 



(2.7) 



r 7 -l/2 m k i m 2 

Z $,k —f ^ 4 ^°"0,fc 

&0,k 



7 In thermal equilibrium at high temperature this decoupling is not obvious. We will comment on this 
point in section 6. 

8 We note that the expressions (2.7) obey the well known Gcll-Mann-Oakes-Renner relation m^f^ = 
-2m(^V) +0{rh 2 ) [33]. 
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Here we have defined the dimensionless, renormalized couplings 

A fc = Z$ 2 k (p = 2al k ) , 



dp 2 vr u ' fc/ ' (2.8) 
hk = Z$J Z^ k h k . 

We will mainly be interested in the "physical values" of the quantities (2.7) in the limit 
k — > where the infrared cutoff is removed, i.e. f w = fn,k=o, m 1 = m/^ k=0 , etc. We point 
out that the formalism of composite fields provides the link [26] to the chiral condensate 
(ip^ since the expectation value a is related to the expectation value of a composite 
quark-ant iquark operator. 



2.2 Flow equations and infrared stability 

At first sight, a reliable computation of T k ^ seems a very difficult task. Without a 
truncation T k is described by an infinite number of parameters (couplings, wave function 
renormalizations, etc.) as can be seen if Tk is expanded in powers of fields and deriva- 
tives. For instance, the sigma mass is obtained as a zero of the exact inverse propagator, 
lim fc ^ T k 2 \q) | $=($), which formally receives contributions from terms in T k with arbitrarily 
high powers of derivatives and the expectation value a . Realistic nonperturbative trunca- 
tions of T fc which reduce the problem to a manageable size are crucial. We will follow here 
a twofold strategy: 

• Physical observables like meson masses, decay constants, etc., can be expanded in 
powers of (current) quark masses in a similar way as in chiral perturbation theory [34]. 
To a given finite order of this expansion only a finite number of terms of a simultaneous 
expansion of Tk in powers of derivatives and $ are required if the expansion point is 
chosen properly. Details of this procedure and some results can be found in [35, 36, 37]. 

• Because of an approximate partial IR fixed point behavior of the flow equations in the 
symmetric regime, i.e. for k x ss < k < k$, the values of many parameters of T k for 
k — > will be almost independent of their initial values at the compositeness scale k$. 
For large enough h k<6 only a few relevant parameters need to be computed accurately 
from QCD. They can alternatively be determined from phenomenology. Because of 
the present lack of an explicit QCD computation we will pursue the latter approach. 

In combination, these two points open the possibility for a perhaps unexpected degree of 
predictive power within the linear quark meson model. We wish to stress, however, that a 
perturbative treatment of the model at hand, e.g., using perturbative RG techniques, cannot 
be expected to yield reliable results. The renormalized Yukawa coupling is very large at the 
scale kq,. Even the IR value of hk is still relatively big 

h k=0 = — ^ ~ 6.5 (2.9) 

Jtt 

and hk increases with k. The dynamics of the linear quark meson model is therefore clearly 
nonperturbative for all scales k < k$. 
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We will now turn to the flow equations for the linear quark meson model. We first note 
that the flow equations for T k and T k — | / d 4 xtr + &jj are identical. The source term 
therefore does not need to be considered explicitly and only appears in the condition (2.5) 
for ($). It is convenient to work with dimensionless and renormalized variables therefore 
eliminating all explicit /c-dependence. With 

u(t,p) = k~ d U k {p) , p = Z^ k k 2 ' d p (2.10) 

and using (1.6) as a first truncation of the effective average action T k one obtains the flow 
equation (t = \n(k/kq>)) 

d 

—u = —du + (d — 2 + 77$) pu' 

dt ( I , (2-11) 

+ 2v d |3/ V; V*) + lo(u' + 2pu"; 77 $ ) - 2* +1 NJ { F )d (-ph 2 ; % ) J . 

Here v^ 1 = 2 d+1 ir d / 2 T(d/2) and primes denote derivatives with respect to p. The number of 
quark colors is denoted as N c . We will always use in the following N c = 3. Eq. (2.11) is a 
partial differential equation for the effective potential u(t,p) which has to be supplemented 
by the flow equation for the Yukawa coupling and expressions for the anomalous dimensions 
77$, r/ip. The symbols l d , l^ d denote bosonic and fermionic mass threshold functions, re- 
spectively, which are defined in appendix B. They describe the decoupling of massive modes 
and provide an important nonperturbative ingredient. For instance, the bosonic threshold 
functions 

l*( w . - n + S n,Q -i h 2n- d f d d q 1 8R k 1 

involve the inverse average propagator P(q 2 ) = q 2 + Z^ l k R k (q 2 ) where the infrared cutoff is 
manifest. These functions decrease ~ w ~( n+1 ) for w 1. Since typically w = M 2 /k 2 with 
M a mass of the model, the main effect of the threshold functions is to cut off fluctuations of 
particles with masses M 2 3> k 2 . Once the scale k is changed below a certain mass threshold, 
the corresponding particle no longer contributes to the evolution and decouples smoothly. 

The dimensionless renormalized expectation value k ee 2k 2 ~ d Zq> tk a 2 )k) with Wo tk the k- 
dependent VEV of $, may be computed for each k directly from the condition (2.5) 

u'(t, re) = -JLzzAT^Zsj/ 2 ee e g . (2.13) 



'2k 

Note that re = in the symmetric regime for vanishing source term. Equation (2.13) allows 
us to follow the flow of re according to 

j t K = 7^{^- d - 2 ^- 2 m u ' M } (2 - 14) 



with A ee u"(t, re). We define the Yukawa coupling for p = re and its flow equation reads [27] 

(2.15) 



— h 2 = (d - 4 + 2r]^, + 1]$) h 2 - 2v d h 4 <3l[ F 1 B)d ( y ^h 2 K, e 9 ; rfy, 77$) 



l[ F i B)d (^h 2 K, e g + 2Are; 77^, 77$) 



15 



Similarly, the the scalar and quark anomalous dimensions are infered from 



+ 2^N c h 2 m[ F)d (h 2 K;r ] ^ , 
7fy = -^\nZ^ :k = 2^h 2 l3m[ F 2 B)d (^h 2 K,e g -,7] i ,,7]^) 



(2.16) 



+ m 



i F 2 B)d {^h 2 K,e g + 2Ak; 77^,77$) | , 



which is a linear set of equations for the anomalous dimensions. The threshold functions 

4f,n2 d > m l,n 2 ) m'p d an d m m^2 are a l so specified in appendix B. 

The flow equations (2.11), (2.14) — (2.16), constitute a coupled system of ordinary and 
partial differential equations which can be integrated numerically. Here we take the effective 
current quark mass dependence of h^, Zq>^ and Z^^ into account by stopping the evolution 
according to eqs. (2.15), (2.16), evaluated for the chiral limit, below the pion mass m n . (For 
details of the algorithm used here see refs. [38, 39].) One finds for d — 4 that chiral symmetry 
breaking indeed occurs for a wide range of initial values of the parameters including the 
presumably realistic case of large renormalized Yukawa coupling and a bare mass of 
order 100 MeV. Driven by the strong Yukawa coupling, the renormalized mass term u'(t, p = 
0) decreases rapidly and goes through zero at a scale k x sB not far below k$. Here the 
system enters the spontaneously broken regime and the effective average potential develops 
an absolute minimum away from the origin. The evolution of of the potential minimum 
a ok = k 2 n/2 turns out to be reasonably stable already before k ~ m n where it stops. 
We take this result as an indication that our truncation of the effective action leads at 
least qualitatively to a satisfactory description of chiral symmetry breaking. The reason 
for the relative stability of the IR behavior of the VEV (and all other couplings) is that 
the quarks acquire a constituent mass M q = ha ~ 300 MeV in the spontaneously broken 
regime. As a consequence they decouple once k becomes smaller than M q and the evolution 
is then dominantly driven by the light Goldstone bosons. This is also important for our 
approximation of neglecting the residual gluonic interactions in the quark sector of the model 
as outlined in subsection 2.1. 

Most importantly, one finds that the system of flow equations exhibits an approximate 
IR fixed point behavior in the symmetric regime [27]. To see this explicitly we study the flow 
equations (2.11), (2.14) — (2.16) subject to the condition (2.2). For the relevant range of p 
both u'(t, p) and u'(t, p) + 2pu"(t, p) are then much larger than ph 2 (t) and we may therefore 
neglect in the flow equations all scalar contributions with threshold functions involving these 



16 



large masses. This yields the simplified equations (d = 4, v 4 1 = 32-7T 2 ) 



d 



-4u + (2 + m )pu'-^J F)4 ( 1 -ph"), 



dt - ■ v- ■ 2vr 2 u v 2 

^ = (2.17) 

% = • 

Of course, it should be clear that this approximation is only valid for the initial range of 
running below k<$> before the (dimensionless) renormalized scalar mass squared u'(t,p = 0) 
approaches zero near the chiral symmetry breaking scale. The system (2.17) is exactly 
soluble. Using m 4 (0) = 1 which holds independently of the choice of the IR cutoff we 
find 

K\t) = Zi 1 (t) = 



h 2 



1 (2.18) 
u(t,p) = e -V(e*^) - 7^ jldre-H^\\h\t)pe*) . 

(The integration over r on the right hand side of the solution for u can be carried out by 
first exchanging it with the one over momentum implicit in the definition of the threshold 
function Zq^ 4 (see appendix B).) Here ui(p) = u(0, p) denotes the effective average potential 
at the compositeness scale and tij is the initial value of h 2 at i.e. for t — 0. For simplicity 
we will use an expansion of the initial value effective potential ui(p) in powers of p around 
p = 

MP) = E ^~P n (2-19) 

even though this is not essential for the forthcoming reasoning. Expanding also Z^ 4 in 
eq. (2.18) in powers of its argument one finds for n > 2 

!gM _ e 2(^ 2)t «W0) iV c (-ir(n-l)! (F)4 r _ 2(n _ 2)t] 
W*(f) " e /zf + n 2 2^ 2 (n-2) ln [[) n L 6 J' [Z - M) 

For decreasing t — > — oo the initial values become rapidly unimportant and u^/h 2n 
approaches a fixed point. For n = 2, i.e., for the quartic coupling, one finds 

^M-.-^Sl (221) 



leading to a fixed point value (u^ /h 2 )* = 1. As a consequence of this fixed point behavior 
the system looses all its "memory" on the initial values u^ 2 ^ at the compositeness scale 
This typically happens before the approximation u'(t,p),u'{t,p) + 2pu"{t,p) ^> ph 2 (t) 
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breaks down and the solution (2.18) becomes invalid. Furthermore, the attraction to partial 
infrared fixed points continues also for the range of k where the scalar fluctuations cannot 
be neglected anymore. The initial value for the bare dimensionless mass parameter 

<(0) _ m k (9 99) 

h 2 ~ P ^ > 

lb j hjQ 

is never negligible. (In fact, using the values for m| # and kq> computed in ref. [26] one 
obtains rn 2 k<6 /k% ~ 0.036.) For large hi (and dropping the constant piece «/(0)) the solution 
(2.18) therefore behaves with growing \t\ as 

Z*(t) ~ -^t, 

, 2M 8tt 2 , 

h\t) ~ -— , (2.23) 

u(t,p) - ^e- 2 V(^-^/^re- 4 ^ 4 (^ 2 Wpe 2 0. 

In other words, for /if — > oo the IR behavior of the linear quark meson model will depend 
(in addition to the value of the compositeness scale kq> and the quark mass rh) only on one 
parameter, m 2 ,^. We have numerically verified this feature by starting with different values 

for u) (0). Indeed, the differences in the physical observables were found to be small. This 
IR stability of the flow equations leads to a perhaps surprising degree of predictive power! 
For definiteness we will perform our numerical analysis of the full system of flow equations 
(2.11), (2.14) — (2.16) with the idealized initial value ui(p) = u'j(0)p in the limit h 2 - — > oo. 
It should be stressed, though, that deviations from this idealization will lead only to small 
numerical deviations in the IR behavior of the linear quark meson model as long as the 
condition (2.2) holds, say for hj <: 15 [27]. 

With this knowledge at hand we may now fix the remaining three parameters of our 
model, m 2 ^ and rh by using f n = 92.4 MeV, the pion mass M n = 135 MeV and the 
constituent quark mass M q as phenomeno logical input. Because of the uncertainty regarding 
the precise value of M q we give in table 2 the results for several values of M q . The first line of 
table 2 corresponds to the choice of M q and A/ = u"(k) which we will use for the forthcoming 
analysis of the model at finite temperature. As argued analytically above the dependence on 
the value of A/ is weak for large enough hj as demonstrated numerically by the second line. 
Moreover, we notice that our results, and in particular the value of j, are rather insensitive 
with respect to the precise value of M q . It is remarkable that the values for k$ and 
are not very different from those computed in ref. [26]. As compared to the analysis of 
ref. [27] the present truncation of T k is of a higher level of accuracy: We now consider 
an arbitrary form of the effective average potential instead of a polynomial approximation 
and we have included the pieces in the threshold functions which are proportional to the 
anomalous dimensions. It is encouraging that the results are rather robust with respect to 
these improvements of the truncation. 

Once the parameters k$, and rh are fixed there are a number of "predictions" of 
the linear meson model which can be compared with the results obtained by other methods 
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M q 
MeV 


A/ 


MeV 


■> 

m, 

hi 
K <s> 


,1/3 

MeV 


m(fc<j) 
McV 


m(l GeV) 
MeV 


(i/>V>)(lGeV) 
McV 3 


r(0) 
Jit 

MeV 


303 


1 


618 


0.0265 


66.8 


14.7 


11.4 


-(186) 3 


80.8 


300 





602 


0.026 


66.8 


15.8 


12.0 


-(183) 3 


80.2 


310 





585 


0.025 


66.1 


16.9 


12.5 


-(180) 3 


80.5 


339 





552 


0.0225 


64.4 


19.5 


13.7 


-(174) 3 


81.4 



Table 2: The table shows the dependence on the constituent quark mass M q of the input parameters 
m^ 4 /fc| and j as well as some of our "predictions". The phenomenological input used here besides M q is 
/tt = 92.4 MeV, = 135 MeV. The first line corresponds to the values for M q and A/ used in the remainder 
of this work. The other three lines demonstrate the insensitivity of our results with respect to the precise 
values of these parameters. 



or direct experimental observation. First of all one may compute the value of m at a scale 
of 1 GeV which is suitable for comparison with results obtained from chiral perturbation 
theory [40] and sum rules [41]. For this purpose one has to account for the running of this 
quantity with the normalization scale from kq> as given in table 2 to the commonly used value 
of 1 GeV: 777.(1 GeV) = A' 1 rri(kq,). A reasonable estimate of the factor A is obtained from 
the three loop running of rh in the MS scheme [41] . For M q ~ 300 MeV corresponding to 
the first two lines in table 2 its value is A ~ 1.3. The results for rh(l GeV) are in acceptable 
agreement with recent results from other methods [40, 41] even though they tend to be 
somewhat larger. Closely related to this is the value of the chiral condensate 



(1 GeV) = -Ami [U Z l]i=o - 2m] • (2.24) 



These results are quite non-trivial since not only f w and ttT| enter but also the computed 
IR value Z<s> t k=o- We emphasize in this context that there may be substantial corrections 
both in the extrapolation from k§ to 1 GeV and in the factor a q (see (A.l)). The latter is 
due to the neglected influence of the strange quark which may be important near These 
uncertainties have only little effect on the physics at lower scales as relevant for our analysis 
of the temperature effects. Only the value of j which is fixed by m w enters here. 

A further more qualitative test concerns the mass of the sigma resonance or radial mode 
whose renormalized mass squared may be computed according to (2.7) in the limit k — > 0. 
From our numerical analysis we obtain Xk=o — 20 which translates into m a ~ 430 MeV. One 
should note, though, that this result is presumably not very accurate as we have employed 
in this work the approximation of using the Goldstone boson wave function renormalization 
constant also for the radial mode. Furthermore, the explicit chiral symmetry breaking 
contribution to is certainly underestimated as long as the strange quark is neglected. 
In any case, we observe that the sigma meson is significantly heavier than the pions. This 
is a crucial consistency check for the linear quark meson model. A low sigma mass would 
be in conflict with the numerous successes of chiral perturbation theory [34] which requires 
the decoupling of all modes other than the Goldstone bosons in the IR-limit of QCD. The 
decoupling of the sigma meson is, of course, equivalent to the limit A — > oo which formally 
describes the transition from the linear to the non-linear sigma model and which appears 
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to be reasonably well realized by the large IR-values of A obtained in our analysis. We also 
note that the issue of the sigma mass is closely connected to the value of fj®, the value 
of f n in the chiral limit rh = also given in table 2. To lowest order in (f n — ffl)/ or, 
equivalently, in rh one has 

f" - = JTjh = ^ ■ < 2 - 25 > 

A larger value of m a would therefore reduce the difference between ffl and /„-. 

In figure 4 we show the dependence of the pion mass and decay constant on the average 
current quark mass rh. These curves depend very little on the values of the initial parameters 
as demonstrated in table 2 by f^°\ We observe a relatively large difference of 12MeV 
between the pion decay constants at rh — rh phys and rh — 0. According to (2.25) this 
difference is related to the mass of the sigma particle and will be modified in the three 
flavor case. We will later find that the critical temperature T c for the second order phase 
transition in the chiral limit is almost independent of the initial conditions. The values of 
and T c essentially determine the non-universal amplitudes in the critical scaling region 
(cf. section 5). In summary, we find that the behavior of our model for small k is quite 
robust as far as uncertainties in the initial conditions at the scale k$ are concerned. We will 
see that the difference of observables between non-vanishing and vanishing temperature is 
entirely determined by the flow of couplings in the range < k ^ 3T. 
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3 Thermal equilibrium and dimensional reduction 



The extension of flow equations to thermal equilibrium situations at non-vanishing temper- 
ature T is straightforward [42]. In the Euclidean formalism non-zero temperature results 
in (anti-)periodic boundary conditions for (fermionic) bosonic fields in the Euclidean time 
direction with periodicity 1/T [43]. This leads to the replacement 



in the trace in (1.1) when represented as a momentum integration, with a discrete spectrum 
for the zero component 

(i\ — j 2/7rT for bosons , , 

qo{L) ~\ (21 + 1)ttT for fermions . {6 } 

Hence, for T > a four-dimensional QFT can be interpreted as a three-dimensional model 
with each bosonic or fermionic degree of freedom now coming in an infinite number of copies 
labeled by I e (Matsubara modes). Each mode acquires an additional temperature 
dependent effective mass term q\{$). In a high temperature situation where all massive 
Matsubara modes decouple from the dynamics of the system one therefore expects to observe 
an effective three-dimensional theory with the bosonic zero modes as the only relevant 
degrees of freedom. In other words, if the characteristic length scale associated with the 
physical system is much larger than the inverse temperature the compactified Euclidean 
"time" dimension cannot be resolved anymore. This phenomenon is known as "dimensional 
reduction" [44]. 

The formalism of the effective average action automatically provides the tools for a 
smooth decoupling of the massive Matsubara modes as the scale k is lowered from k 3> T 
to fc < T. It therefore allows us to directly link the low-T, four-dimensional QFT to 
the effective three-dimensional high-T theory. The replacement (3.1) in (1.1) manifests 
itself in the flow equations (2.11), (2.14) — (2.16) only through a change to T-dependent 
threshold functions. For instance, the dimensionless functions l^iw] 77$) defined in eq. (2.12) 
are replaced by 

l d (w T - t*) = n + S <* v -W-*r V f ( 1 dRkiq2) ) 1 ( 3 3) 



1&2Z- 



where q 2 = q% + q 2 and q = 2nlT. A list of the various temperature dependent threshold 
functions appearing in the flow equations can be found in appendix C. There we also 
discuss some subtleties regarding the definition of the Yukawa coupling and the anomalous 
dimensions for T 7^ 0. In the limit k ^> T the sum over Matsubara modes approaches the 
integration over a continuous range of q and we recover the zero temperature threshold 
function l^(w;r]^). In the opposite limit k <^ T the massive Matsubara modes (/ 7^ 0) are 
suppressed and we expect to find a d — 1 dimensional behavior of l„. In fact, one obtains 
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from (3.3) 

l d n (w,T/k;^) ~ l d n {w- m ) for T < fc , 

l d (w,T/k;r}$) ~ Z^V;^) for T > fc . 

For our choice of the infrared cutoff function eq. (1.3), the temperature dependent 
Matsubara modes in l d (w,T/k;r]^) are exponentially suppressed for T k whereas the 
behavior is more complicated for other threshold functions appearing in the flow equations 
(2.11), (2.14) — (2.16). Nevertheless, all bosonic threshold functions are proportional to 
T/k for T 3> k whereas those with fermionic contributions vanish in this limit 9 . This 
behavior is demonstrated in figure 5 where we have plotted the quotients lf(w,T/k)/lf(w) 
and l[ F) \w,T/k)/l{ F) \w) of bosonic and fermionic threshold functions, respectively. One 
observes that for k ^> T both threshold functions essentially behave as for zero temperature. 
For growing T or decreasing k this changes as more and more Matsubara modes decouple 
until finally all massive modes are suppressed. The bosonic threshold function if shows for 
k <^ T the linear dependence on T/k derived in eq. (3.4). In particular, for the bosonic 
excitations the threshold function for w C 1 can be approximated with reasonable accuracy 
by ln(w; 77$) for T/k < 0.25 and by (4T / k)l^(w; 77$) for T/k > 0.25. The fermionic threshold 
function Z^ 4 tends to zero for k C T since there is no massless fermionic zero mode, i.e. in 
this limit all fermionic contributions to the flow equations are suppressed. On the other 
hand, the fermions remain quantitatively relevant up to T/k ~ 0.6 because of the relatively 
long tail in figure 5b. The transition from four to three-dimensional threshold functions 
leads to a smooth dimensional reduction as k is lowered from k ^> T to k <C T! Whereas for 
k ^> T the model is most efficiently described in terms of standard four-dimensional fields 
$ a choice of rescaled three-dimensional variables $3 = $/\/T becomes better adapted for 
k <T. Accordingly, for high temperatures one will use a potential 

k k 
u 3 (t,p 3 ) = -u(t,p) ; p 3 = —p . (3.5) 

In this regime r fe ^ corresponds to the free energy of classical statistics and r fc>0 is a classical 
coarse grained free energy. 

For our numerical calculations at non-vanishing temperature we exploit the discussed 
behavior of the threshold functions by using the zero temperature flow equations in the range 
k > 10T. For smaller values of k we approximate the infinite Matsubara sums (cf. eq. (3.3)) 
by a finite series such that the numerical uncertainty at k — 10T is better than 10~ 4 . This 
approximation becomes exact in the limit k <C 10T. 



4 The quark meson model at T ^ 

In section 2.1 we have considered the relevant fluctuations that contribute to the flow of 
in dependence on the scale k. In a thermal equilibrium situation Y k also depends on the 

9 For the present choice of Rk the temperature dependence of the threshold functions is considerably 
smoother than in ref. [42]. 
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Figure 5: The plot shows the temperature dependence of the bosonic (a) and the fermionic (b) threshold 
functions lf(w,T/k) and l[ F ^ 4 (w,T /fc), respectively, for different values of the dimensionless mass term w. 
The solid line corresponds to w = whereas the dotted ones correspond to w = 0.1, w = 1 and w = 10 
with decreasing size of the dots. For r>l the bosonic threshold function becomes proportional to T/k 
whereas the fermionic one tends to zero. In this range the theory with properly rescaled variables behaves 
as a classical three-dimensional theory. 
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temperature T and one may ask for the relevance of thermal fluctuations at a given scale 
k. In particular, for not too high values of T (cf. sect. 6) the "initial condition" for the 
solution of the flow equations should essentially be independent of temperature. This will 
allow us to fix T k<s> from phenomenological input at T = and to compute the temperature 
dependent quantities in the infrared (k — > 0). We note that the thermal fluctuations which 
contribute to the r.h.s. of the flow equation for the meson potential (2.11) are effectively 
suppressed for T ^ k/A (cf. section 3). Clearly for T <; k^>/3 temperature effects become 
important at the compositeness scale. We expect the linear quark meson model with a 
compositeness scale k<$> ~ 600 MeV to be a valid description for two flavor QCD below a 
temperature of about 10 170 MeV. 

We compute the quantities of interest for temperatures T ^ 170 MeV by solving nu- 
merically the T-dependent version of the flow equations (2.11), (2.14) — (2.16) (cf. section 3 
and appendix C) by lowering k from k§ to zero. For this range of temperatures we use 
the initial values as given in the first line of table 2. This corresponds to choosing the zero 
temperature pion mass and the pion decay constant (f n = 92.4 MeV for m n = 135 MeV) as 
phenomenological input. The only further input is the constituent quark mass M q which 
we vary in the range M q ~ 300 — 350 MeV. We observe only a minor dependence of our 
results on M q for the considered range of values. In particular, the value for the critical 
temperature T c of the model remains almost unaffected by this variation. 

We have plotted in figure 6 the renormalized expectation value 2<x of the scalar field 
as a function of temperature for three different values of the average light current quark 
mass rh. (We remind that 2<7o(T = 0) = f^.) For rh = the order parameter do of chiral 
symmetry breaking continuously goes to zero for T — > T c = 100.7 MeV characterizing the 
phase transition to be of second order. The universal behavior of the model for small T — T c 
and small rh is discussed in more detail in section 5. We point out that the value of T c 
corresponds to ffi = 80.8 MeV, i.e. the value of the pion decay constant for rh — 0, which 
is significantly lower than f n = 92.4 MeV obtained for the realistic value rn p hy S . If we would 
fix the value of the pion decay constant to be 92.4 MeV also in the chiral limit (rh = 0), 
the value for the critical temperature would raise to 115 MeV. The nature of the transition 
changes qualitatively for m^fl where the second order transition is replaced by a smooth 
crossover. The crossover for a realistic m ph ys or m^T = 0) = 135 MeV takes place in a 
temperature range T ~ (120 — 150) MeV. The middle curve in figure 6 corresponds to a 
value of rh which is only a tenth of the physical value, leading to a zero temperature pion 
mass m n = 45 MeV. Here the crossover becomes considerably sharper but there remain 
substantial deviations from the chiral limit even for such small quark masses rh ~ 1 MeV. 
The temperature dependence of m n has already been mentioned in the introduction (see 
fig. 2) for the same three values of rh. As expected, the pions behave like true Goldstone 
bosons for rh — 0, i.e. their mass vanishes for T < T c . Interestingly, m n remains almost 
constant as a function of T for T < T c before it starts to increase monotonically. We 
therefore find for two flavors no indication for a substantial decrease of m n around the 

10 There will be an effective temperature dependence of induced by the fluctuations of other degrees of 
freedom besides the quarks, the pions and the sigma which are taken into account here. We will comment on 
this issue in section 6. For realistic three flavor QCD the thermal kaon fluctuations will become important 
for T k, 170 MeV. 
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Figure 6: The expectation value 2cto is shown as a function of temperature T for three different values of 
the zero temperature pion mass. 



critical temperature. 

The dependence of the mass of the sigma resonance m a on the temperature is displayed 
in figure 7 for the above three values of rh. In the absence of explicit chiral symmetry 
breaking, rh — 0, the sigma mass vanishes for T < T c . For T <T C this is a consequence of 
the presence of massless Goldstone bosons in the Higgs phase which drive the renormalized 
quartic coupling A to zero. In fact, A runs linearly with k for T k/4 and only evolves 
logarithmically for T ^ k/4. Once ??i^0 the pions acquire a mass even in the spontaneously 
broken phase and the evolution of A with k is effectively stopped at k ~ m n . Because of 
the temperature dependence of <7o,fc=o (cf. figure 6) this leads to a monotonically decreasing 
behavior of m a with T for T ^ T c . This changes into the expected monotonic growth once 
the system reaches the symmetric phase for 11 T > T c . For low enough rh one may use the 
minimum of m a (T) for an alternative definition of the (pseudo-) critical temperature denoted 
as Table 1 in the introduction shows our results for the pseudocritical temperature 

for different values of m or, equivalently, m 7r (T = 0). For a zero temperature pion mass 
m n = 135 MeV we find T$ = 128 MeV. At larger pion masses of about 230 MeV we 
observe no longer a characteristic minimum for m a apart from a very broad, slight dip at 
T ~ 90 MeV. A comparison of our results with lattice data is given in section 5. In fig. 8 we 
show the renormalized quartic coupling A as a function of temperature for two fixed values 
of the average current quark mass. The upper curve corresponds to the physical value of rh 
or, equivalently, m^iT = 0) = 135 MeV whereas the lower curve shows A for rh = 0. One 
observes the vanishing of the renormalized quartic coupling in the chiral limit for T < T c as 

11 See section 2.2 for a discussion of the zero temperature sigma mass. 
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Figure 7: The plot shows the m a as a function of temperature T for three different values of the zero 
temperature pion mass. 
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Figure 8: The plot shows the renormalized quartic scalar self coupling A as a function of temperature T 
for the physical value of m (solid line) as well as for m = (dashed line). 
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Figure 9: The plot shows the renormalized $ 6 scalar self coupling U3 as a function of temperature T in 
the chiral limit. 

discussed above. The renormalized scalar <3> 6 self interaction 

U 3 (T) = Z^^^l(p = 2al(T)) (4.1) 

assumes a small negative value for realistic quark masses in the temperature range T ^ 
35 MeV with 2U 3 (T)a$(T) ~ -0.5 < A(T) and 2C/ 3 (T)^(T) ~ 8.0,8.5,1.5 for T = 
80, 120, 160 MeV. We display U 3 (T) in figure 9 for the chiral limit where one observes a 
discontinuity of U 3 (T) at the critical temperature T c . 

Our results for the chiral condensate (ipi/)^ as a function of temperature for different 
values of the average current quark mass are presented in figure 1 in the introduction. We 
will compare (^ip^ (T,m) with its universal scaling form for the 0(4) Heisenberg model in 
section 5. Another interesting quantity is the temperature dependence of the constituent 
quark mass. Figure 10 shows M q (T) for m = 0, rh = m p h ys /10 and m = rh p h ys , respectively. 
Its behavior is related to the temperature dependence of the renormalized order parameter 
<7o(T) = <T fe=0 (T) and the renormalized Yukawa coupling h(T) = hk=o(T). The temperature 
dependence of h in the chiral limit can be found in fig. 11. Near the critical temperature one 
notices a characteristic dip. This results from the long wavelength pion fluctuations through 
a non-analytic behavior of the mesonic wave function renormalization Z$(T) = Z$ t k=o(T) 
which is displayed in figure 12. There we also present the temperature dependence of the 
fermionic wave function renormalization Z^{T) = Z^,^ =0 (T). Away from the chiral limit 
we take the effective quark mass dependence of hk(T), Z$^(T) and Z^^(T) into account 
by stopping their evolution when k reaches the temperature dependent pion mass. In this 
way we observe no substantial quark mass dependence of these quantities except for Z$(T), 
and consequently for h(T), in the vicinity of the critical temperature. A more complete 
truncation would incorporate field dependent wave function renormalization constants and 
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Figure 10: The plot shows the constituent quark mass M q as a function of T for three different values 
of the average light current quark mass rh. The solid line corresponds to the realistic value rh = m p hy S 
whereas the dotted line represents the situation without explicit chiral symmetry breaking, i.e., rh = 0. The 
intermediate, dashed line assumes rh = Wphys/lO. 
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Figure 11: The plot shows the Yukawa coupling, h, as a function of temperature T in the chiral limit. 
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Figure 12: The plot shows the scalar (solid line) and quark (dashed line) wave function renormalization 
constants, Z$(T) and Z^{T) ■ 1CT 2 , respectively, as functions of temperature T for rh = 0. 



a field dependent Yukawa coupling. 

Our ability to compute the complete temperature dependent effective meson potential 
U is demonstrated in fig. 13 where we display the derivative of the potential with respect 
to the renormalized field 4>r = (Z$p/2) 1 / 2 , for different values of T. The curves cover a 
temperature range T = (5— 175) MeV. The first one to the left corresponds to T = 175 MeV 
and neighboring curves differ in temperature by AT = 10 MeV. One observes only a weak 
dependence of dU(T)/d(f)R on the temperature for T ^ 60 MeV. Evaluated for (fi R = <r this 
function connects the renormalized field expectation value with m n (T), the source j and the 
mesonic wave function renormalization Z$(T) according to 



dU{T) 



{<t>R = Co) 



2j 



Z l J 2 (T) 



= Aa (T)ml(T) . 



(4.2) 



We close this section with a short assessment of the validity of our effective quark meson 
model as an effective description of two flavor QCD at non-vanishing temperature. The 
identification of qualitatively different scale intervals which appear in the context of chiral 
symmetry breaking, as presented in section 2.1 for the zero temperature case, can be gener- 
alized to T 7^ 0: For scales below k§ there exists a hybrid description in terms of quarks and 
mesons. For k x sB < k ^ 600 MeV chiral symmetry remains unbroken where the symmetry 
breaking scale k x sB{T) decreases with increasing temperature. Also the constituent quark 
mass decreases with T (cf. figure 10). The running Yukawa coupling depends only mildly 
on temperature for T ^ 120 MeV (see fig. 11). (Only near the critical temperature and for 
rh = the running is extended because of massless pion fluctuations.) On the other hand, 
for k ^ 4T the effective three-dimensional gauge coupling increases faster than at T = 
leading to an increase of Aq CD (T) with T [18]. As k gets closer to the scale Aq CD (T) it is 
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Figure 13: The plot shows the derivative of the meson potential U{T) with respect to the renormalized 
field <j) R = (Z^p/2) 1 / 2 for different values of T. The first curve on the left corresponds to T = 175 MeV. 
The successive curves to the right differ in temperature by AT = 10 MeV down to T = 5 MeV. 

no longer justified to neglect in the quark sector confinement effects which go beyond the 
dynamics of our present quark meson model. Here it is important to note that the quarks 
remain quantitatively relevant for the evolution of the meson degrees of freedom only for 
scales k T/0.6 (cf. fig. 5, section 3). In the limit k <C T/0.6 all fermionic Matsubara modes 
decouple from the evolution of the meson potential according to the temperature dependent 
version of eq. (2.11). Possible sizeable confinement corrections to the meson physics may 
occur if Aqcd(^) becomes larger than the maximum of M q (T) and T/0.6. From fig. 10 
we infer that this is particularly dangerous for small m in a temperature interval around 
T c . Nevertheless, the situation is not dramatically different from the zero temperature case 
since only a relatively small range of k is concerned. We do not expect that the neglected 
QCD non-localities lead to qualitative changes. Quantitative modifications, especially for 
small m and |T — T c \ remain possible. This would only effect the non-universal amplitudes 
(see sect. 5). The size of these corrections depends on the strength of (non-local) deviations 
of the quark propagator and the Yukawa coupling from the values computed in the quark 
meson model. 

5 Universal critical behavior 

In this section we study the linear quark meson model in the vicinity of the critical tem- 
perature T c close to the chiral limit m = 0. In this region we find that the sigma mass m^ 1 
is much larger than the inverse temperature T -1 , and one observes an effectively three- 
dimensional behavior of the high temperature quantum field theory. We also note that the 
fermions are no longer present in the dimensionally reduced system as has been discussed 
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in section 3. We therefore have to deal with a purely bosonic 0(4)-symmetric linear sigma 
model. At the phase transition the correlation length becomes infinite and the effective 
three-dimensional theory is dominated by classical statistical fluctuations. In particular, 
the critical exponents which describe the singular behavior of various quantities near the 
second order phase transition are those of the corresponding classical system. 

Many properties of this system are universal, i.e. they only depend on its symmetry 
(0(4)), the dimensionality of space (three) and its degrees of freedom (four real scalar com- 
ponents). Universality means that the long-range properties of the system do not depend on 
the details of the specific model like its short distance interactions. Nevertheless, important 
properties as the value of the critical temperature are non-universal. We emphasize that 
although we have to deal with an effectively three-dimensional bosonic theory, the non- 
universal properties of the system crucially depend on the details of the four-dimensional 
theory and, in particular, on the fermions. 

Our aim is a computation of the critical equation of state which relates for arbitrary 
T near T c the derivative of the free energy or effective potential U to the average current 
quark mass rh. The equation of state then permits to study the temperature and quark 
mass dependence of properties of the chiral phase transition. 

At the critical temperature and in the chiral limit there is no scale present in the theory. 
In the vicinity of T c and for small enough rh one therefore expects a scaling behavior of the 
effective average potential Uk [45] and accordingly a universal scaling form of the equation 
of state. There are only two independent scales close to the transition point which can be 
related to the deviation from the critical temperature, T — T c , and to the explicit symmetry 
breaking through the quark mass rh. As a consequence, the properly rescaled potential 
can only depend on one scaling variable. A possible choice for the parameterization of the 
rescaled "unrenormalized" potential is the use of the Widom scaling variable [46] 



Here j3 is the critical exponent of the order parameter a in the chiral limit rh = (see 
equation (5.5)). With U'(p = 2al) = j/(2a ) the Widom scaling form of the equation of 
state reads [46] 



where the exponent 5 is related to the behavior of the order parameter according to (5.7). 
The equation of state (5.2) is written for convenience directly in terms of four-dimensional 
quantities. They are related to the corresponding effective variables of the three-dimensional 
theory by appropriate powers of T c . The source j is determined by the average current quark 
mass rh as j = Zm^rh. The mass term at the compositeness scale, ml , also relates the 

chiral condensate to the order parameter according to = ~~ 2m| $ (a — rh) . The critical 

temperature of the linear quark meson model was found in section 4 to be T c = 100.7 MeV. 

The scaling function / is universal up to the model specific normalization of x and itself. 
Accordingly, all models in the same universality class can be related by a rescaling of a 
and T — T c . The non-universal normalizations for the quark meson model discussed here 



x = 




(5.1) 



(2<x /T c ) 




(5.2) 
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are defined according to 

f(0) = D , f(-B-W) = 0. (5.3) 

We find D = 1.82 ■ 10~ 4 , B = 7.41 and our result for (3 is given in table 3. Apart from the 
immediate vicinity of the zero of f(x) we find the following two parameter fit (cf. ref. [47]) 
for the scaling function, 

/ fit (a;) = 1.816 • 1(T 4 (1 + 136.1a;) 2 (1 + 160.9 #:r) A 
(1 + 160.9 (0.9446 0^-i/( 7 -2-a) x y-2-A 

to reproduce the numerical results for / and df/dx at the 1 — 2% level with 6 = 0.625 (0.656), 
A = —0.490 (—0.550) for x > (x < 0) and 7 as given in table 3. The universal properties 
of the scaling function can be compared with results obtained by other methods for the 
three-dimensional 0(4) Heisenberg model. In figure 14 we display our results along with 
those obtained from lattice Monte Carlo simulation [48], second order epsilon expansion 
[49] and mean field theory. We observe a good agreement of average action, lattice and 
epsilon expansion results within a few per cent for T < T c . Above T c the average action and 
the lattice curve go quite close to each other with a substantial deviation from the epsilon 
expansion and mean field scaling function 12 . 

Before we use the scaling function f(x) to discuss the general temperature and quark 
mass dependent case, we consider the limits T = T c and m — 0, respectively. In these 
limits the behavior of the various quantities is determined solely by critical amplitudes and 
exponents. In the spontaneously broken phase (T < T c ) and in the chiral limit we observe 
that the renormalized and unrenormalized order parameters scale according to 

2^ (T) 2(Tc -TY' 2 



(2E)' 



T, 



c 



(5.5) 

2ao(T) b (T c -T^ 



1 



respectively, with E = 0.814 and the value of B given above. In the symmetric phase the 
renormalized mass m = = m a and the unrenormalized mass m = Z^ 2 m behave as 



T r ^ I V T, 



V 

1 



(5.6) 

m(T) _ ^ c+ y^( T - T c\<' 



T r v > \ T, 



where £ + = 0.270, C + = 2.79. For T = T c and non-vanishing current quark mass we have 

(5.7) 



2a u < • n 1/6 




T 

J- c 

with the value of D given above. 

12 We note that the question of a better agreement of the curves for T < T c or T > T c depends on the 
chosen non-universal normalization conditions for x and / (cf. eq. (5.3)). 
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Figure 14: The figure shows a comparison of our results, denoted by "average action" , with results of other 
methods for the scaling function of the three-dimensional 0(4) Heisenberg model. We have labeled the axes 
for convenience in terms of the expectation value ao and the source j of the corresponding four-dimensional 
theory. The constants B and D specify the non-universal amplitudes of the model (cf. eq. 5.3). The curve 
labeled by "MC" represents a fit to lattice Monte Carlo data. The second order epsilon expansion [49] and 
mean field results are denoted by "e" and "mf" , respectively. Apart from our results the curves are taken 
from [48]. 
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5 


P 


V 


average action 


0.787 


1.548 


4.80 


0.407 


0.0344 


e 


0.73(2) 


1.44(4) 


4.82(5) 


0.38(1) 


0.03(1) 


MC 


0.7479(90) 


1.477(18) 


4.851(22) 


0.3836(46) 


0.0254(38) 



Table 3: The table shows the critical exponents corresponding to the three-dimensional (9(4)-Heisenberg 
model. Our results are denoted by "average action" whereas "e" labels the exponents obtained from the 
4 — e expansion to seven loops. The bottom line contains lattice Monte Carlo results. [50] 



Though the five amplitudes E, B, £ + , C + and D are not universal there are ratios of 
amplitudes which are invariant under a rescaling of a and T — T c . Our results for the 
universal amplitude ratios are 

R x = C + DB 5 - 1 = 1.02 , 

R ( = = 0>852 ? (5. 8 ) 

i + E = 0.220 . 

Those for the critical exponents are given in table 3. Here the value of rj is obtained from the 
temperature dependent version of eq. (2.16) (cf. appendix C) at the critical temperature. For 
comparison table 3 also gives the results of the 4— e expansion to seven loops [51, 52] 13 as well 
as lattice Monte Carlo results [50] which have been used for the lattice form of the scaling 
function in fig. 14. 14 There are only two independent amplitudes and critical exponents, 
respectively. They are related by the usual scaling relations of the three-dimensional scalar 
0(iV)-model [52] which we have explicitly verified by the independent calculation of our 
exponents. 

We turn to the discussion of the scaling behavior of the chiral condensate (^pip^ for the 
general case of a temperature and quark mass dependence. In figure 1 in the introduction we 
have displayed our results for the scaling equation of state in terms of the chiral condensate 15 

(Irp3 \ V<5 
fit)) +3 (5 ' 9) 

as a function of T/T c = 1 + x(j/T^f(x)) 1 ^ 5 for different quark masses or, equivalently, 
different values of j. The curves shown in figure 1 correspond to quark masses m = 0, m = 
m p h ys /10, m = m p h y s and m = 3.5m phys or, equivalently, to zero temperature pion masses 
m n = 0, m n = 45MeV, = 135 MeV and m n = 230 MeV, respectively (cf. figure 4). One 
observes that the second order phase transition with a vanishing order parameter at T c for 
m = is turned into a smooth crossover in the presence of non-zero quark masses. 

13 The values are quoted e.g. in ref. [50]. 

14 See also ref. [53] and references therein for a recent calculation of critical exponents using similar methods 
as in this work. For high precision estimates of the critical exponents see also ref. [54, 55]. 

15 In the literature also a different definition of the chiral condensate is used, corresponding to (ipi)) = 
-n4W(T c 3 /(z))] 1/5 . 
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The scaling form (5.9) for the chiral condensate is exact only in the limit T — > T c , 
j — > 0. It is interesting to find the range of temperatures and quark masses for which (iptp^ 
approximately shows the scaling behavior (5.9). This can be infered from a comparison (see 
fig. 1) with our full non-universal solution for the T and j dependence of (ipi^ as described 
in section 4. For m n = one observes approximate scaling behavior for temperatures 
T <: 90 MeV. This situation persists up to a pion mass of m n = 45 MeV. Even for the 
realistic case, = 135 MeV, and to a somewhat lesser extent for m n = 230 MeV the 
scaling curve reasonably reflects the physical behavior for T ^ T c . For temperatures below 
T c , however, the zero temperature mass scales become important and the scaling arguments 
leading to universality break down. 

The above comparison may help to shed some light on the use of universality arguments 
away from the critical temperature and the chiral limit. One observes that for temperatures 
above T c the scaling assumption leads to quantitatively reasonable results even for a pion 
mass almost twice as large as the physical value. This in turn has been used for two 
flavor lattice QCD as theoretical input to guide extrapolation of results to light current 
quark masses. From simulations based on a range of pion masses 0.3 ^ m n /m p ^ 0.7 and 
temperatures < T < 250 MeV a "pseudocritical temperature" of approximately 140 MeV 
with a weak quark mass dependence is reported [56]. Here the "pseudocritical temperature" 
T pc is defined as the inflection point of (ipip^ as a function of temperature. The values of 
the lattice action parameters used in [56] with N t — 6 were am = 0.0125, Q/g 2 = 5.415 and 
am = 0.025, Q/g 2 = 5.445. For comparison with lattice data we have displayed in figure 1 
the temperature dependence of the chiral condensate for a pion mass m n = 230 MeV. From 
the free energy of the linear quark meson model we obtain in this case a pseudocritical 
temperature of about 150 MeV in reasonable agreement with the results of ref . [56] . In 
contrast, for the critical temperature in the chiral limit we obtain T c = 100.7 MeV. This 
value is considerably smaller than the lattice results of about (140 — 150) MeV obtained 
by extrapolating to zero quark mass in ref. [56]. We point out that for pion masses as 
large as 230 MeV the condensate (ipij^ (T) is almost linear around the inflection point for 
quite a large range of temperature. This makes a precise determination of T c somewhat 
difficult. Furthermore, figure 1 shows that the scaling form of (^pip^ (T) underestimates 
the slope of the physical curve. Used as a fit with T c as a parameter this can lead to an 
overestimate of the pseudocritical temperature in the chiral limit. We also mention here the 
results of ref. [57]. There two values of the pseudocritical temperature, T pc = 150(9) MeV 
and T pc = 140(8), corresponding to am = 0.0125, Q/g 2 = 5.54(2) and am = 0.00625, 
Q/g 2 = 5.49(2), respectively, (both for N t = 8) were computed. These values show a 
somewhat stronger quark mass dependence of T pc and were used for a linear extrapolation 
to the chiral limit yielding T c = 128(9) MeV. 

The linear quark meson model exhibits a second order phase transition for two quark 
flavors in the chiral limit. As a consequence the model predicts a scaling behavior near the 
critical temperature and the chiral limit which can, in principle, be tested in lattice simu- 
lations. For the quark masses used in the present lattice studies the order and universality 
class of the transition in two flavor QCD remain a partially open question. Though there 
are results from the lattice giving support for critical scaling [58, 59] there are also recent 
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simulations with two flavors that reveal significant finite size effects and problems with 0(4) 
scaling [60, 61]. 

6 Additional degrees of freedom 

So far we have investigated the chiral phase transition of QCD as described by the linear 
0(4)-model containing the three pions and the sigma resonance as well as the up and 
down quarks as degrees of freedom. Of course, it is clear that the spectrum of QCD is 
much richer than the states incorporated in our model. It is therefore important to ask to 
what extent the neglected degrees of freedom like the strange quark, strange (pseudo) scalar 
mesons, (axial) vector mesons, baryons, etc., might be important for the chiral dynamics 
of QCD. Before doing so it is perhaps instructive to first look into the opposite direction 
and investigate the difference between the linear quark meson model described here and 
chiral perturbation theory based on the non-linear sigma model [34]. In some sense, chiral 
perturbation theory is the minimal model of chiral symmetry breaking containing only the 
Goldstone degrees of freedom. By construction it is therefore only valid in the spontaneously 
broken phase and can not be expected to yield realistic results for temperatures close to T c 
or for the symmetric phase. However, for small temperatures (and momentum scales) the 
non-linear model is expected to describe the low-energy and low-temperature limit of QCD 
reliably as an expansion in powers of the light quark masses. For vanishing temperature it 
has been demonstrated recently [36, 37] that the results of chiral perturbation theory can 
be reproduced within the linear meson model once certain higher dimensional operators in 
its effective action are taken into account for the three flavor case. Moreover, some of the 
parameters of chiral perturbation theory (L 4 , . . . , L 8 ) can be expressed and therefore also 
numerically computed in terms of those of the linear model. For non-vanishing temperature 
one expects agreement only for low T whereas deviations from chiral perturbation theory 
should become large close to T c . Yet, even for T <T C small quantitative deviations should 
exist because of the contributions of (constituent) quark and sigma meson fluctuations in 
the linear model which are not taken into account in chiral perturbation theory. 

From [62] we infer the three-loop result for the temperature dependence of the chiral 
condensate in the chiral limit for N light flavors 

{N 2 — IT 2 N 2 — 1 / T 2 \ 2 

+ ^-''(iS?) ln F7} + 0(T8 »- 

The scale Ti can be determined from the D-wave isospin zero tttt scattering length and 
is given by Ti = (470 ± 100) MeV. The constant F is (in the chiral limit) identical to 
the pion decay constant F = /£>) = 80.8 MeV. In fi gure 15 we have plotted the chiral 
condensate as a function of T/Fq for both, chiral perturbation theory according to (6.1) 
and for the linear quark meson model. As expected the agreement for small T is very 
good. Nevertheless, the anticipated small numerical deviations present even for T <^ T c 
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Figure 15: The plot displays the chiral condensate (ipip) as a function of T/fk°\ The solid line corre- 
sponds to our results for vanishing average current quark mass m = whereas the dashed line shows the 
corresponding three-loop chiral perturbation theory result for Ti = 470 MeV. 

due to quark and sigma meson loop contributions are manifest. For larger values of T, 
say for T ^ 0.8/^ the deviations become significant because of the intrinsic inability of 
chiral perturbation theory to correctly reproduce the critical behavior of the system near 
its second order phase transition. 

Within the language of chiral perturbation theory the neglected effects of thermal quark 
fluctuations may be described by an effective temperature dependence of the parameter 
F (T). We notice that the temperature at which these corrections become important equals 
approximately one third of the constituent quark mass M q (T) or the sigma mass m a (T), 
respectively, in perfect agreement with fig. 5. As suggested by this figure the onset of the 
effects from thermal fluctuations of heavy particles with a T-dependent mass m H {T) is 
rather sudden for T <: m,H(T)/3. These considerations also apply to our two flavor quark 
meson model. Within full QCD we expect temperature dependent initial values at kq>. 

The dominant contribution to the temperature dependence of the initial values presum- 
ably arises from the influence of the mesons containing strange quarks as well as the strange 
quark itself. Here the quantity m| seems to be the most important one. (The temperature 
dependence of higher couplings like A(T) is not very relevant if the IR attractive behavior 
remains valid, i.e. if remains small for the range of temperatures considered. We ne- 

glect a possible T-dependence of the current quark mass m.) In particular, for three flavors 
the potential contains a term 

(det $ + det $ f ) = -V kit ip s <$> uu <$> dd + ... (6.2) 
which reflects the axial C/a(1) anomaly. It yields a contribution to the effective mass term 
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proportional to the expectation value ($ ss ) = (fi s , i-e. 

A ™L = ~\ v k^s ■ (6.3) 

Both, X7 fc<t and y9 s , depend on T. We expect these corrections to become relevant only for 
temperatures exceeding m K (T)/3 or M 8 (T)/3. We note that the temperature dependent 
kaon and strange quark masses, m^(T) and M S (T), respectively, may be somewhat different 
from their zero temperature values but we do not expect them to be much smaller. A 
typical value for these scales is around 500 MeV. Correspondingly, the thermal fluctuations 
neglected in our model should become important for T ^ 170 MeV. It is even conceivable 
that a discontinuity appears in y? s (T) for sufficiently high T (say T ~ 170 MeV). This 
would be reflected by a discontinuity in the initial values of the 0(4)-model leading to a 
first order transition within this model. Obviously, these questions should be addressed 
in the framework of the three flavor SUl{3) x SUr(3) quark meson model. Work in this 
direction is in progress. 

We note that the temperature dependence of v(T)ip s (T) is closely related to the question 
of an effective high temperature restoration of the axial symmetry [3, 8]. The rj 

mass term is directly proportional to this combination, m^(T) — m^(T) ~ ^v(T)tp s (T) 
[35]. Approximate Ua(1) restoration would occur if <p s (T) or v(T) would decrease sizeable 
for large T. For realistic QCD this question should be addressed by a three flavor study. 
Within two flavor QCD the combination u k (p s is replaced by an effective anomalous mass 
term . The temperature dependence of (T) could be studied by introducing quarks 
and the axial anomaly in the two flavor matrix model of ref. [39]. We add that this question 
has also been studied within full two flavor QCD in lattice simulations [60, 63, 64]. So far 
there does not seem to be much evidence for a restoration of the Ua(X) symmetry near T c 
but no final conclusion can be drawn yet. 

To summarize, we have found that the effective two flavor quark meson model presum- 
ably gives a good description of the temperature effects in two flavor QCD for a temperature 
range T ^ 170 MeV. Its reliability should be best for low temperature where our results 
agree with chiral perturbation theory. However, the range of validity is considerably ex- 
tended as compared to chiral perturbation theory and includes, in particular, the critical 
temperature of the second order phase transition in the chiral limit. We have explicitly 
connected the universal critical behavior for small |T — T c \ and small current quark masses 
with the renormalized couplings at T = and realistic quark masses. The main quantita- 
tive uncertainties from neglected fluctuations presumably concern the values of and T c 
which, in turn, influence the non-universal amplitudes B and D in the critical region. We 
believe that our overall picture is rather solid. Where applicable our results compare well 
with numerical simulations of full two flavor QCD. 
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Appendices 



A The quark mass term 

In this appendix we determine the source j = diag(j u , jj, . . .) as a function of the average 
current quark mass rh. In this context it is important to note that the source j does not 
depend on the IR cutoff scale k. Since j is determined by the properties of the quark 
meson model at the compositeness scale kq> and also enters directly the value of the pion 
mass, which is determined at k — 0, this relation provides a bridge between the short and 
long distance properties of the quark meson model. This will allow us to compute the 
chiral condensate (vpip^ or the parameter Bq of chiral perturbation theory [34]. (We expect, 
however, sizeable corrections when going from two to three flavors. They arise because of 
the relevance of strange quark physics at scales near In a more general context we 

need the proportionality coefficient a q between the source j q and the current quark mass 
m q , q — u,d, . . ., taken at the renormalization scale 16 /x = 

Jq = ^a q m q . (A.l) 

For a computation of the coefficient a q we need to look into the details of the introduction 
of composite meson fields in QCD [26, 28]. Let us assume that at the scale fc$ a part of the 
QCD average action for quarks r fc4 ,[?/>] factorizes in the quark bilinear 

j^y9{P,<l)^Lb(p)^Ra{p + q) (A.2) 

such that 

r k M = -F k M + nM- ( A -3) 

We can then introduce meson fields by inserting the identity 

TV J £>$exp(-F fc , + = 1 (A.4) 

into the path integral which formally defines r fc<i) [?/>]. (Here AMs a field independent nor- 
malization factor.) This effectively replaces in (A. 3) the term —F k<ji [x] by 17 

W + + *] = gM*.fe) + l ax Ztiw) ^ mM) + ' ' ' ' (A - 5) 



16 We will occasionally use the notation m(it), m q (fi) or (tptp) (/z) = ('0V , ) fe= o (m) ( not to be confused with 
(tptp). = (ip^u (M = k<s>)) in order to indicate the renormalization scale fi. If no argument is given \i = fc$ 
is assumed. 

17 Thc summation over internal indices as well as the integration over momenta has been suppressed. For 
complex Xab(q) similar terms have to be supplemented in the expansion. See ref. [26, 28] for a more detailed 
description. In our Euclidean conventions one has Xat ~ +9*i } m^La- 
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The original multi-quark interaction —F k<s> [x] is canceled by the lowest order term in the 
Taylor expansion in $. Instead, we have substituted mesonic self-interactions i*fc # [3>] and 
interactions between mesons and quarks corresponding to the terms in the expansion which 
contain powers of x an d In particular, we may specialize to the case where the derivative 
terms in F k<s> are small and consider a local form F k<s> = J d 4 xf k9 (x)- A quark mass term is 
linear in x an d translates into a source term for $ 



2 5 '2 

where m = diag(m u ,m rf , . . .). A factorizing four fermion interaction yields 



k Tr X h — ^,Tr$t$ + m^Tr^x + x t$) . ( A J) 

The second term corresponds to the Yukawa interaction with h k , = 7n k<s g- We can therefore 
extract a q from eq. (A. 6) as 

« ? = H- (A-8) 
We note that only the terms linear and quadratic in x influence the value of a q . We could 
either restrict the composite fields from the beginning to the ones contained in the 0(4)- 
symmetric linear a-model or work with all the fields contained in a complex 2x2 matrix $. 
In the latter case the anomaly term would contribute to both the masses and the Yukawa 
coupling. The net result is the same with m\ denoting the relevant mass term for the 0(4) 
vector. For our conventions with h k<s> = 1 we have to normalize with g = rn k ~ 2 . Finally an 
eight fermion interaction becomes 

^X k , Tr( X t x )2 _^ I A fe jTr$t$) 2 + A fc ,Tr$t$ Tr( X t$ + $t % ) + . . . . (A . 9) 

We see here the appearance of terms quadratic in the quarks involving higher powers of <3>. 

There is an alternative, equivalent way of understanding the relation between j and m q : 
The quark masses in the picture with mesons must be equal at the scale k§ to the current 
quark mass m q (kq>). Let us consider an 0(4)-symmetric fermionic interaction m| Tr^x + 
|Ajfc # (Trx^x) 2 which leads to a meson potential 

U k ,=ml^T&<$>+ l -\ k , (Tr$t$) 2 . (A.IO) 

In the mesonic picture the quarks acquire masses through the Yukawa coupling to $ 

h 



p-(i + ^{*%{*) k ){*) k (A.n; 
Am \ m k ) 



where the second term arises from the higher order coupling in (A. 9). Here ($) fc = 
diag(<y9 M , ifd-, . . .) is the expectation value at the coarse graining scale k in the presence of the 
source term and M k = diag(M u ,Mrf, . . .). It is sufficient to specify the dependence of U k , 
on real diagonal fields Then the ip q are determined from the condition 

g^M = 2 {ml, + A fe4) £ \ { fq=3q- (A.12) 
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Identifying M k=k<Sf in (A. 11) with m{k§) one has 

««(i + |r£«) =^ = 2< + 2^Ef4 ( A - 13 ) 

and we recover (A. 8) or, in our normalization with = 1, h k<s , = 1, 

j = 2m^m. (A. 14) 

It is remarkable that higher order terms (e.g. ~ A^) do not influence the relation between 
j and m. Only the quadratic term m| enters which is in our scenario the only relevant 
coupling. This feature is an important ingredient for the predictive power of the model as 
far as the absolute size of the current quark mass is concerned. An appearance of higher 
order couplings in a q would make it very hard to compute this quantity. We emphasize that 
within our formalism there is no difference of principle between the current quark mass and 
the constituent quark mass. Whereas the current quark mass m q (k<$>) at the normalization 
scale /i — k$ corresponds to M q;k at the compositeness scale /c$ the constituent quark mass 
is M q ^=o- As k is lowered from k$ to zero one observes that the quark mass increases, 
similarly to the running current quark mass. Once chiral symmetry breaking sets in at the 
scale k x sB there is a large increase in the quark masses, especially for M u and 

The formalism of composite fields also provides the link [26] to the chiral condensate 
since the expectation value (<&) is related to the expectation value of a composite 
quark- ant iquark operator. For A = one has [28] 

+ ($t) fc = __L (m) (A;*) + m,(fc*) + m\{k*) (A.15) 
with (ipip) a suitably regularized operator normalized at ji — 



B Threshold functions 

In this appendix we list the various definitions of dimensionless threshold functions appear- 
ing in the flow equations and the expressions for the anomalous dimensions for T = 0. They 
involve the inverse scalar average propagator P(q) defined in (1.4) and the corresponding 
fermionic function P F which can be chosen as [27] 

P F ( q ) = P( q )= q 2 (l + r F (q)) 2 . (B.l) 

We abbreviate 

1 = ' 2 • p ^ s p w • ^ ik p{x) ■ wt p s rJ t p • (R2) 

etc., and use the formal definition 

d_ J_dRkd_ 2 P F djZ^rp] d 

dt ~ Z^ ik dt dP Z^ k 1 + r F dt dP F ' 1 ' ) 
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(B.4) 



The bosonic threshold functions read 

#(iu;7fr) = l d (w) - r)J d n {w) 

2 Jo \Z^ k dt J V J 

l d n u n 2 {Wl,W 2 ]m) = lZ una (w 1 ,W2)-rfri* ltn2 (w 1 ,W 2 ) 

= _^2(n 1+ n 2 )-^ ° dx J-l^_^p + Wlh ?y ni (p + ^ 

where n,rii,n 2 > is assumed. For n^O the functions may also be written as 
The fermionic integrals l^ d (w; rf^) = l^p d (w) — rj^l^ d (w) are defined analogously as 

fi.») = (n + s„, )e-* r^x*- i ~ p : 8lz * a : rF] ( p +^r" >+ " • 

jo i + 7> or x 7 

(B.6) 

Furthermore one has 

CSVi, ^2; rfc) = i^K, w 2 ) - ^2) - WK>i, ^2) 

2 Jo 9t \ 



P F (x) + + £; 2 w 2 ]' 

"ini^^i.^;^*) = m^ ljn2 («;i,«;2) - ri*m* ltna (w 1 ,w 2 ) 

JO 



2 Jo 9t ( [P(x) + k 2 wi] ni [P(x) + fc 2 w 2 ] n2 J 

mi (w; = m\ [wj - r]^m\ [wj 

2 Jo 9t V<9xP F (x) + A; 2 w; 

^2; Tfy, = mJ^K, «> 2 ) - ^™SVi, «> 2 ) - fte^XVi, w 2 ) 

= _I^2(n 1+ n 2 -l)- d r dxx jll ] 

2 Jo 9t \ [Pf(x) + A; 2 wi]"i [P(x) + A; 2 w 2 ]™2 J ' 

(B.7) 

The dependence of the threshold functions on the anomalous dimensions arises from the 
t-derivative acting on and Z^^ within R k and Z^ ;k r F , respectively. We furthermore 
use the abbreviations 

i d M = i d n {^m) , 4 F)d (%)^4 F)d (0;%) 
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etc. and note that in four dimensions the integrals 

Z 2 4 (0, 0) = if )4 (0, 0) = Zi™(0, 0) = mf )4 (0) = m£ B)4 (0, 0) = 1 (B.9) 
are independent of the particular choice of the infrared cutoff. 

C Temperature dependent threshold functions 

Non-vanishing temperature manifests itself in the flow equations (2.11), (2.14) — (2.16) only 
through a change to T-dependent threshold functions. In this appendix we will define these 
functions and discuss some subtleties regarding the definition of the anomalous dimensions 
and the Yukawa coupling for T ^ 0. The corresponding T = threshold functions can be 
found in appendix B where also some of our notation is fixed. 

The flow equation (2.11) for the effective average potential involves a bosonic and a 
fermionic threshold function whose generalization to finite temperature is straightforward 

J d (wf-n.) - + ^ n ' o) Vd ~ l k 2n - d+1 f V I'" dr r^Z- 1 dtRM 

l^ d (w,f; V *) = (n + M— (C.l) 

x f E f^^ZTi Pf(vf) d t [Z,, k r F (y F )} 
f^h *> k l + r F (y F ) [P F (y F )+wk 2 ] n+1 



(C.2) 



where f = T/k and 

y = x + (2/ttT) 2 

y F = x + (2/ + l)VT 2 

The computation of the anomalous dimensions 77$, 77^ and the flow equation for the 
Yukawa coupling h at non-vanishing temperature requires some care. The anomalous di- 
mensions determine the IR cutoff scale dependence of Z®^ and Z^^ according to 77$ = 
— dt In Z^, 77^, = —dt In with t = Ink/ k$>. It is important to realize that for a computa- 
tion of the scale dependence of the effective three-dimensional Zq>^ and Z^ for T 7^ mo- 
mentum dependent wave function renormalization constants of the four-dimensional theory 
are required. This is a consequence of the fact that in the three-dimensional model each of 
the infinite number of different Matsubara modes of a four-dimensional bosonic or fermionic 
field (f>(Q) corresponds to a different value of Q = 2nlT or Q = (21 + l)vrT, respectively, 
with Q 2 = Ql + Q 2 and I G ZZ . We will therefore allow for momentum dependent wave 
function renormalizations, i.e. for a kinetic part of of the form 

^ = / { Z *> kiq2)q2Tl O^fo)) + ZtAfMQWQMQ)} ( c -3) 
in momentum space. 
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In the 0(4)-model the evolution equation for Z<$>^{Q) m &y then be obtained by consid- 
ering a background field configuration with a small momentum dependence, 



$j(x) = (pSj! + (5^e- iQx + 8<p*e iQx ) 8 j2 ; j = 1, . . . , 4 



(C.4) 
supplemented by 

i; a =^ a = 0; a=l,2. (C.5) 

Expanding around this configuration at the minimum of the effective average potential Uk 
we observe that Sip corresponds to an excitation in the Goldstone boson direction. The 

(2) 

exact inverse two-point function rjj. turns out to be block-diagonal with respect to scalar 
and fermion indices for this configuration. It therefore decays into corresponding matrices 



T^l and acting in the scalar and fermion subspaces, respectively. The scale dependence 
of the scalar wave function renormalization for non-vanishing T is obtained from (1.1) and 
(C.3) for the configuration (C.4) as 



,(2) 



d_ 
dt 



z$,k{Q 2 



l 



lim 



8 



- Tr 



Q2 ys<ps<p*^o 8(8ip8(p' 
i d 



Tr 



(2) 
Sk 



R, 



-id_ 

dt 



Rk 



(Tpl + R Fk ) 



dt 



R 



Fk 



(C.6) 



(Q->0) 



In the three-dimensional theory there is now a different scalar wave function renormalization 
Z$,k(Qo,Q) f° r each Matsubara mode Qq. As in the four-dimensional model for T = 
we neglect the momentum dependence of the wave function renormalization constants 
and evaluate Z$>^ for Q = for each Matsubara mode. We will furthermore simplify 
the truncation of the effective average action by choosing the Matsubara zero-mode wave 
function renormalization constant for all Matsubara modes, i.e., approximate 



Z^ k (f) = Z^k(Ql = 0,Q 2 = 0) . 



(C.7) 



This is justified by the rapid decoupling of all massive Matsubara modes within a small 
range of k for fixed T as discussed in section 3. This results in the expression (2.16) for 77$ 
but now with temperature dependent threshold functions (T = T/k) 



_ 1 7 2(m+n 2 -l)-d+l dl! d~l 

~ 2 



(d-l)v d 



Hv) 



Hv) 



doc cc ^ \ 

dt \ [P(y) + k 2 wi] n ^ [P(y) + k 2 w 2 ] n2 



AF)d 



5-d 



dv d -i 



2" (d-l)v d le2Z 
For further technical details we refer the reader to ref. [27]. 



did l + r F (y F ) 

axx 2 y F - 



dt \dx P F (y F ) + k 2 w / 



(C.8) 
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The fermion anomalous dimension and the flow equation for the Yukawa coupling can 
be obtained by considering a field configuration 



$j(x) = (pSj!-, j = l,...,4 

Mx) = ^e- tQx (C.9) 
i, a {x) = ^ a e iQx ; a = 1,2. 

The derivation follows similar lines as for the scalar anomalous dimension discussed above. 
For computational details we refer the reader to [27]. An important difference as compared 
to Z<p t k(Q) relates to the fact that there are no fermionic zero modes. It would therefore 
be inconsistent to define Z^iT) or hk(T) at Qo = if Q denotes the external fermion 
momentum. Yet, we will again resort to the approximation of using the same wave function 
renormalization constant and Yukawa coupling for all fermionic Matsubara modes. For 
the same reason as for Z^ t k{T) we will use for this purpose the mode with the lowest In- 
dependent mass, i.e. define 

Z^, k (f) = Z^ k (Ql = ir 2 T 2 ,Q 2 = 0) 
h(f) = h k (Ql = n 2 T 2 ,Q 2 = 0), 

where we have neglected a possible dependence of h k on the external scalar momentum of 
the Yukawa vertex. This yields the expressions (2.15) and (2.16) for the flow of h 2 and rj^, 
respectively, but now with the T-dependent threshold functions 



rn 



[ F 2 B)d (w u W2,T; 71^,7]$) = m[ F 2 B)d (w 1 ,w 2 ,T) 



- V^i F 2 B)d {wi,w 2 , T) - %ra^ S)d (u;i, w 2 , T) 

2 (d-l)v d ^Jo 

x gf l + r F (y F ) P(y) \ 

dt \ \P F {yF) + Pwi]™ 1 [P(y) + k 2 w 2 ] n 'i J 

= K2{n 1+ n 2 )- d+ lVj-l f>< y, [°° dx x <¥ 

2 v d l£Z ,Jo 

II I 1. 

dt \ [Pf(j/f) + k 2 wi} n ^[P{y) + k 2 w 2 ] n * J 



(C.11) 
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